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Abstract 

We present ground-based and HST optical and infrared observations of Swift XRF 1003 16D / SN 2010bh. It 
is seen that the optical light curves of SN 2010bh evolve at a faster rate than the archetype GRB-SN 1998bw, 
but at a similar rate to SN 2006aj, a supernova that was spectroscopically linked with XRF 060218, and at a 
similar rate to non-GRB associated type Ic SN 19941. We estimate the rest-frame extinction of this event from 
our optical data to be E(B -V) = 0.18 ± 0.08 mag. We find the V-band absolute magnitude of SN 2010bh to 
be My = — 18.62 ± 0.08, which is the faintest peak V-band magnitude observed to-date for a spectroscopically- 
confirmed GRB-SNe. When we investigate the origin of the flux at t — 1 = 0.598 days, it is shown that the 
light is not synchrotron in origin, but is likely coming from the supernova shock break-out. We then use 
our optical and infrared data to create a quasi-bolometric light curve of SN 2010bh which we model with 
a simple analytical formula. The results of our modeling imply that SN 2010bh synthesized a nickel mass 
of Mn; ss O.1M , ejected M e j ~ 2.2M and has an explosion energy of E^ m 1.4 x 10 52 erg. Thus, while 
SN 2010bh is an energetic explosion, the amount of nickel created during the explosion is much less than that 
of SN 1998bw, and only marginally more than SN 19941. Finally, for a sample 22 GRB-SNe we check for a 
correlation between the stretch factors and luminosity factors in the R band and conclude that no statistically- 
significant correlation exists. 

Subject headings: Gamma Rays: Bursts — Stars: Supernovae: General — Stars: Supernovae: Individual: 
Alphanumeric: SN 2010bh 



1. INTRODUCTION 

For more than a decade the intriguing connection between 
long-duration gamma ray bursts (GRBs) and X-ray Flashes 
(XRFs) with type Ibc, stripped-envelope, core-collapse super- 
novae (SNe) has been studied in considerable depth. What 
is now the archetype gamma ray burst-supernova (GRB-SN), 
SN 1998bw (e.g. Galama et al. 1998a; Patat et al. 2001) 
provided the first observational evidence that GRBs are gen- 
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erated during the collapse of a massive star. Since then sev- 
eral other detections of GRB-SNe have strengthened the link 
between the collapse of massive stars and GRBs, including: 
GRB 030329 / SN 2003dh (Hjorth et al. 2003; Stanek et al. 
2003; Matheson et al. 2003), GRB 031203 / SN 20031w (e.g. 
Malesani et al. 2004; Gal- Yam et al. 2004) and XRF 060218 
/ SN 2006aj (e.g. Pian et al. 2006; Campana et al. 2006; 
Ferrero et al. 2006; Mazzali et al. 2006). In all of the afore- 
mentioned cases, spectroscopy of each event revealed unusu- 
ally broad spectral features, indicative of much higher ejecta 
velocities that are much higher than those seen in "typical" Ic 
SN events. 

However, there exists not just a simple dichotomy between 
typical type Ic SNe such as SN 19941 (e.g. Richmond et 
al. 1996; Iwamoto et al. 1994) that have modest peak 
magnitudes and ejecta moving at velocities anywhere from 
8,000- 18,000 km s~ ! (e.g. Matheson et al. 2001), and 
GRB-SNe that generally have brighter peak magnitudes and 
ejecta velocities up to 30,000 km s -1 (e.g. SN 1998bw; 
Galama et al. 1998a). Indeed the type Ic SN class is quite 
heterogeneous, especially when the discoveries of broad-lined 
SNe that are not accompanied by a GRB trigger are also con- 
sidered. Events such as SN 2007 gr (Paragi et al. 2010; though 
see Soderberg et al. 2010a) and SN 2009bb (Stritzinger et al. 
2009; Pignata et al. 2011) reveal that relativistic ejecta are 
present in some non-GRB associated SNe. Even more intrigu- 
ingly, Soderberg et al. (2010b) showed that SN 2009bb most 
likely had a prolonged central engine. The multi-frequency 
radio observations presented by Soderberg et al. (2010b) 
are well described by a self-absorbed synchrotron spectrum 
which they suggest is produced as the shock break-out accel- 
erated electrons in the local circumstellar medium. 
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Thus, it seems that the special requisites for GRB- 
production are very atypical of the general Ic SNe population; 
indeed it is possible to have an energetic type Ic SN with a 
bright maximum, broad spectral features, and prolonged cen- 
tral engine activity, without the production of a GRB. Not 
suprisingly, the differences between GRB-SNe and local Ic 
SNe without a GRB-trigger remains elusive, though the in- 
depth investigation of individual cases provides an excellent 
opportunity for understanding these enigmatic events. Previ- 
ous investigations into individual events has revealed that no 
two GRB-SNe are the same, and the latest example we present 
in this paper is no exception. 

XRF 100316D was detected by Swift on the 16 th March, 
2010 at 12:44:50 UT. The redshift was measured to be z = 
0.0591 ±0.0001 (Vergani et al. 2010a; Starling et al. 2011; 
Chornock et al. 2010), and we use this value throughout the 
paper. Spectroscopy presented by Chornock et al. (2010) 
revealed similar behaviour to previous GRB-SNe, where the 
spectra are featureless during the first few days, and with the 
emergence of broad spectral features after several days. As for 
other GRB-SNe, Chornock et al. (2010) reported that there 
was no evidence for helium in their spectroscopy, indicating 
that the progenitor is a highly-stripped star. However, spec- 
troscopy taken by Bufano et al. (2011) revealed a broad He 
absorption feature at 10830A, a spectral line that has only be 
observed one other time in a GRB-SNe (SN 1998bw, Patat 
et al. 2001). Bufano et al. (201 1) also detected He in their 
early spectra at 5876A. These detections of He indicate that 
the progenitor star still retained at least a small fraction of its 
He envelope prior to explosion, indicating that SN 2010bh is 
a type Ibc SNe. 

Starling et al. (2011) found that the high-energy proper- 
ties of XRF 100316D / SN 2010bh are quite different from 
those of other GRB-SNe, but are similar to the high energy 
properties of XRF 060218 / SN 2006aj. Starling et al. (201 1) 
showed in the X-ray spectrum at 144 s < t — 1 < 737 s the 
presence of a soft, hot, black-body component that contributes 
w 3% to the total 0.3 - 10.0 keV flux (though see also Fan et 
al. 2011 who claim that both a Cutoff Power-law model or 
a Cutoff Power-law model + black-body component provide 
equally acceptable fits). The presence of the black-body com- 
ponent (kT = 0.14 keV; T « 1.6 x 10 6 K) is similar in tem- 
perature to the soft, hot, black-body component seen in the 
X-ray spectrum of XRF 060218 / SN 2006aj (kT w 0.17 keV; 
r w 2.0 x 10 6 K; Campana et al. 2006). 

A detailed discussion of the host galaxy metallicity at the 
explosion site of SN 2010bh was performed by Levesque et 
al. (2011). The authors presented evidence, which was also 
reported by Chornock et al. (2010) and Starling et al. (201 1), 
that SN 2010bh occurred in a star-forming region of the host 
galaxy that has a star-formation rate of ~ 1.7M yr _1 (at 
the site of the supernova) and low metallicity (Z < 0.4Z fi ; 
Chornock et al. (2010), Levesque et al. (201 1)). 

In this paper we present in Section [2] optical and in- 
frared photometry of XRF 100316D / SN 2010bh obtained 
on Faulkes Telescope South (FTS); Gemini South (Gemini- 
S) and the Hubble Space Telescope (HST). In Section [3] we 
use our optical light curves (LCs) to derive the peak optical 
properties of SN 2010bh, showing that SN 2010bh evolves 
much quicker than SN 1998bw, and at a rate similar to 
SN 2006aj. We also discuss in Section 3 the origin of the flux 
at t — 1 = 0.598 days, concluding that it is not synchrotron 
but is likely coming from the shock break-out. In Section |4] 



we make an estimate of the rest frame extinction, and derive 
peak absolute magnitudes of SN 2010bh. In Section [5] we 
use our optical and infrared (IR) photometry, along with our 
estimation of the host extinction, to construct the bolometric 
light curve, and use an analytical model to extract physical 
parameters of the SN. In Section|6]we discuss SN 2010bh in 
relation to SN 2006aj, as well as other GRB-SNe and local 
type Ic SNe that were not accompanied by a GRB. In the dis- 
cussion we investigate the light curve properties of most of the 
spectroscopically and photometrically-linked GRB-SNe. By 
making simple assumptions, we obtain host- and afterglow- 
subtracted "supernova" light curves, and compare them with 
SN 1998bw, deriving stretch (s) and luminosity (k) factors of 
the GRB-SNe in relation to SN 1998bw. When then check 
for a correlation between the stretch and luminosity factors, 
concluding that no statistically-significant correlation exists. 

Throughout the paper observer-frame times are used un- 
less specified otherwise in the text. Foreground reddening of 
E(B — V) = 0.117 mag has been corrected for using the dust 
maps of Schlegel, Finkbeiner & Davis (1998). We adopt a flat 
ACDM cosmology with H = 71 km s" 1 Mpc" 1 , i2 M = 0.27, 
and £2a = 1 — D.m = 0.73. For this cosmology a redshift of 
z = 0.0591 corresponds to a luminosity distance of dL =261 
Mpc and a distance modulus of 37.08 mag. 

2. OBSERVATIONS & PHOTOMETRY 
2.1. Data Acquisition & Calibration 
Faulkes Telescope South & Gemini-South 

We obtained forty-one epochs of photometric data between 
March and May 2010 using the 2m Faulkes Telescope South 
(FTS). Observations were made using the 4.6' x 4.6' field of 
view Merope Camera in filters BVR c i. A 42'"' and final epoch 
of observations were made in December 2010 in all filters, 
which were used as template/reference images for image sub- 
traction (see Section 2.2). Images were also taken of Landolt 
photometric standard regions (Landolt 1992) in BVR c i on the 
2 nd of May, 2010 using the Spectral Camera, which with its 
wider field of view (10.5' x 10.5'), allowed us to calibrate 
over 140 secondary standards in the field of XRF 1003 16D. 
Subsequent observations made on Gemini-South and HST are 
calibrated to these secondary standards. 

Five epochs of photometric data were obtained on Gemini- 
South H (Gemini-S) in filters griz, the first starting only rs 
0.5d after the initial burst. A final epoch was obtained on 
the 28''' of January, 2011 in all filters, which were used as 
template images for image subtraction (see Section 2.2). 

Hubble Space Telescope 

We obtained five epochs of photometric data with HST, 
using the Wide Field Camera 3 (WFC3). The first three 
epochs of data were obtained in filters F336W, F555W, 
F8UW, F125W and F160W, while the fourth epoch yielded 
images in filters F336W, F555W, F8UW, and the fifth 
epoch in filters F555VK and F814W. We find for the opti- 
cal transient (OT) associated with XRF 100316D a position 
of: 07 /, 10'"30.54' s (±0.02) s , -56 rf 15'l9.80"(±0.10)". 

19 The data presented here supersedes those published in the GCN Circu- 
lars: GCNs 10513 (Vergani et al. 2010b); 10523 (Levan et al. 2010); 10525 
(Wiersema et al. 2010). 
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FIG. 1.— UBVR C UH light curves of XRF 100316D /SN 2010bh (solid points (ground-based data) and filled squares (HST)). Plotted for comparison are the 
light curves of two known GRB-SNe: SN 1998bw (Galama et al. 1998a; Patat et al. 2001 ; Sollerman et al. 2000) and SN 2006aj (Campana et al. 2006; Sollerman 
et al. 2006; Ferrero et al. 2006), as well as a local SN that is not associated with a GRB: SN 19941 (Richmond et al. 1996). All of the SN 2010bh light curves 
have been corrected for foreground and host extinction. The light curves of the comparison SNe have been shifted by different amounts in each filter to match 
in peak brightness. No other changes have been made to the light curves apart from the light curves of SN 1998bw which have been stretched in each filter by a 
factor of ~ 0.6 (see Table|5} and the early-time data (t — 1„ < 3 days) for SN 2006aj, for which the light was dominated by flux coming from the shock break-out, 
was not included to avoid confusion between different filters. 



Aperture photometry was performed on the HST images us- 
ing standard routines in IRAFPT. A small aperture was used, 
and an aperture correction for an isolated star was computed 
and applied. The aperture-corrected, instrumental magni- 
tudes were then calibrated via standard star photometry to sec- 
ondary standards in the field of XRF 1003 16D. Images taken 
in filters F555W and F814W are calibrated respectively to 
filters V and i using a zero-point and a colour term. The re- 
maining filters are calibrated to AB magnitudes using the ap- 
propriate AB zero-point for each filter that are listed in the 
WFC3 instrument and data-handboolQ The HST magni- 
tudes, which have been corrected for foreground extinction, 
are listed in Table [TJ The quoted errors, which are added in 
quadrature, are derived from the uncertainties associated with 
the photometry and calibration. 

2.2. Image Subtraction 

20 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 

21 http://www.stsci.edu/hst/wfc3 



Using ISIS (Alard 2000), we performed image subtraction 
on the FTS and Gemini-S images, using the last epoch of im- 
ages taken in each filter and on each telescope as the respec- 
tive templates. Subtracting the template from the early FTS 
images gave a clear detection of the optical transient (OT) and 
an accurate position, and the same was also observed for the 
Gemini-S images. 

Aperture photometry was then performed on the subtracted 
images using a small aperture and an aperture correction was 
computed from the pre-subtracted images and applied. The 
instrumental magnitudes obtained from the FTS subtracted 
images are calibrated to BVR c i. The Gemini-S subtracted im- 
ages taken in filters griz are calibrated to BR c iz respectively. 
Images in filters gr are calibrated using a zero-point and a 
colour term, while images taken in i are calibrated using only 
a zero-point. Magnitudes of the secondary standards in z are 
computed using transformation equations from Jester et al. 
(2005), and the Gemini-S data are calibrated to these using a 
zero-point. All of the image-subtracted, foreground-corrected 
magnitudes obtained from the FTS and Gemini-S data are dis- 
played in Table |2j where the quoted errors (added in quadra- 
ture) are determined from the photometry and calibration. 
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TABLE 1 

HST Photometry of XRF 1003 16D / SN 2010bh 



Filter 


t — tr. ( davs) 


Mag 


(7 (maff) 

w filing/ 


Calibrated to 


F336W 


8.650 


21.86 


0.06 


AB 


F336W 


15.810 


22.80 


0.08 


AB 


F336W 


24.860 


23.56 


0.11 


AB 


F336W 


47.340 


24.05 


0.14 


AB 


y 


8.663 


19.50 


0.01 


Vega 


\> 


15.822 


19.81 


0.01 


Vega 


\r 


24.869 


20.52 


0.02 


Vega 


\> 


47.352 


21.63 


0.03 


Vega 


w 
\ 


1 37 490 


23 07 


0.08 


Vega 




8.675 


19.00 


0.01 


AB 




15.844 


19.12 


0.01 


AB 




24.881 


19.59 


0.01 


AB 




47.370 


20.87 


0.02 


AB 


( 


137.518 


22.70 


0.06 


AB 


F125W 


8.720 


19.89 


0.07 


AB 


FU5W 


15.870 


20.09 


0.07 


AB 


F125W 


25.480 


20.10 


0.07 


AB 


F160W 


8.730 


19.93 


0.08 


AB 


F160W 


15.880 


20.18 


0.08 


AB 


F160W 


25.520 


20.32 


0.07 


AB 



All magnitudes have been corrected for foreground extinction. 
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FIG. 2.— Colour curves of XRF 100316D/SN 2010bh, which are corrected 
for foreground and host extinction. Plotted for comparison are the colour 
curves of known GRB-SNe: 1998bw and 2006aj as well as local type Ic 
SN 19941 that was not accompanied by a GRB trigger. No other changes 
have been made to the colour curves apart from those of SN 1998bw, which 
have been stretched by s = 0.63 in B — V and V — R and by s = 0.6 1 in R c — i 
(see Section 3.5). The stretch factors, which are the average values of s deter- 
mined in the individual filters, have been applied so as to present a consistent 
analysis of the properties between SN 2010bh and SN 1998bw. 



2.3. Host Photometry 

As already noted by Starling et al (2011), the host galaxy 
system of XRF 100316D / SN 2010bh has a highly disturbed 
morphology, with a bright central region, a possible spiral 
arm and a number of bright knots (see Figure 8 in Starling 
et al. 2011). SN 2010bh is located on top of a bright knot 
close to the bright central region, which is perhaps the galac- 
tic nucleus. In our final epoch of Gemini-S images, the host- 
complex appears as a blended, elliptical object due to the 
lower resolution of Gemini-S relative to HST. 

We have performed photometry on the final Gemini-S im- 
ages using SExtractor (Bertin & Arnouts 1996). We obtained 
(1) corrected isophotal magnitudes (MAG JSOCOR), and (2) 
Kron-like elliptical aperture magnitudes (MAG_AUTO), with 
the best of these values being used (MAG_BEST). The host 
magnitudes in filters griz have been calibrated to secondary 
standards in the field of XRF 1003 16D, where the magni- 
tudes of the secondary standards in filters grz are computed 
using transformation equations from Jester et al. (2005), and 
the host photometry are calibrated to these using a zero-point. 
The Gemini-S photometry, which have been corrected for 
foreground extinction, are listed in Table [3] where the quoted 
errors of the Gemini-S magnitudes are added in quadrature 
and are derived from the photometry and calibration. 

3. OPTICAL & INFRARED LIGHT CURVES AND 
COLOUR CURVES 

3.1. Light Curves 

The image-subtracted light curves of XRF 100316D / 
SN2010bh, which have been corrected for foreground and 
host extinction (see Section 4), are displayed in Figure[T] Plot- 
ted for comparison are the multi-wavelength light curves of 
two known GRB-SNe: SN 1998bw and SN 2006aj, as well 
as SN 19941 which is a local type Ic SN that was not asso- 
ciated with a GRB. All of the comparison light curves have 
been shifted in magnitude to match SN 2010bh around max- 
imum light. No other changes have been made to the light 



curves apart from the light curves of SN 1998bw, which have 
been stretched in each filter by the values listed in Table [5] 
(see Section 3.5), and the early-time data (t — t D < 3 days) for 
SN 2006aj, for which the light was dominated by flux com- 
ing from the very bright shock break-out, was not included to 
avoid confusion between different filters. The stretch-factors 
have been applied to SN 1998bw for clarity. 

The peak times and magnitudes in each filter are estimated 
from low-order polynomial fits. We have determined that 
SN 2010bh peaked at t-t = 8.02±0.12 days in the B fil- 
ter, and at subsequently later times in the redder filters. The 
quoted errors are statistical and are determined from the dis- 
tribution of peak times determined from the different order 
polynomials. The behaviour of the LC peaking at later times 
in redder filters is typical behaviour for a core-collapse SNe, 
and SN 2010bh appears no different in this respect. The main 
photometric parameters of XRF 100316D / SN 2010bh in fil- 
ters BVR c i are presented in Table |4] 

Comparing the optical light curves of SN 2010bh with those 
of SN 2006aj, it appears they evolve at roughly the same rate 
in the B filter. However it appears that the f/-band light curve 
of SN 2006aj evolves faster than SN 2010bh, while in the 
redder V, R c and i filters it evolves more slowly. It is also 
worth noting that the f/-band measurement of SN 2010bh at 
t — to = 24.8 days is w 0.5 mag brighter than the shifted LC 
of SN 2006aj, and perhaps even more at t — 1 =47.3 days. 
As the fluxes are measured in a fixed, and albeit small, aper- 
ture, the flux we are detecting is from both the SN and the 
underlying host galaxy. When performing photometry on the 
third and fourth epoch [/-band images, it appears that the 
background/host flux is becoming dominant and affecting the 
overall shape of the LC, (i.e. making it appear that it is fading 
slower than is expected). 

3.2. Colour Curves 

The colour curves of SN 2010bh are plotted in Figure 
|2j and have been corrected for foreground and host extinc- 
tion (see Section 4). Plotted for comparison are the colour 
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TABLE 2 

Ground-based Photometry of XRF 1003 16D / SN 2010bh 



Filter 


t — (days) 


Mag 


<y (mag) 


Calibrated to 


Telescope 17 


Filter 


t — ty (days) 


.X- 1 . 1 i: 


a (mag) 


Calibrated to 


Telescope 


B 


0.598 


20.45 


0.03 


Vega 


GS 


R c 


1 1.015 


19.33 


0.02 


Vega 


FTS 


B 


3.489 


20.53 


0.01 


Vega 


GS 


R t 


1 1.985 


19.32 


0.04 


Vega 


FTS 


B 


5.467 


20.22 


0.01 


Vega 


GS 


Rc 


16.045 


19.49 


0.04 


Vega 


FTS 


B 


7.061 


20.25 


0.05 


Vega 


FTS 


Rc 


17.023 


19.52 


0.02 


Vega 


FTS 


B 


7.940 


20. 14 


0.03 


Vega 


FTS 


R t 


19.922 


19.71 


0.02 


Vega 


FTS 


B 


9.956 


20. 1 3 


0.03 


Vega 


FTS 


R c 


">") CV7Q 

LL» /8 


19.83 


0.02 


Vega 


FTS 


B 


10.965 


20.19 


0.06 


Vega 


FTS 


^, 


25.956 


20.07 


0.03 


Vega 


FTS 


B 


16.987 


21.18 


0.07 


Vega 


FTS 


R c 


29.941 


20.54 


0.04 


Vega 


FTS 


B 


17.863 


21.30 


0.05 


Vega 


FTS 


R< 


33.901 


20.83 


0.04 


Vega 


FTS 


B 


19.958 


21.44 


0.13 


Vega 


FTS 


R c 


36.901 


21.08 


0.07 


Vega 


FTS 


B 


22.994 


21.83 


0. 1 1 


Vega 


FTS 




37.988 


21.20 


0. 1 3 


Vega 


FTS 


B 


25.955 


21.88 


0.16 


Vega 


FTS 


Rc 


45.921 


21.36 


0.08 


Vega 


FTS 


B 


30.974 


22.27 


0.10 


Vega 


FTS 


Rc 


48.867 


21.45 


0.07 


Vega 


FTS 


B 


33.910 


22.09 


0.14 


Vega 


FTS 


R c 


54.853 


21 .57 


0.07 


Vega 


FTS 


V 


5.007 


19.64 


0.02 


Vega 


FTS 


R c 


57.358 


21.64 


0.07 


Vega 


FTS 


V 


7.085 


19.49 


0.03 


Vega 


FTS 


1 


0.479 


20.04 


0.01 


AB 


GS 


V 


7.956 


19.48 


0.02 


Vega 


FTS 


1 


3.497 


19.44 


0.01 


AB 


GS 


V 


8.937 


19.47 


0.02 


Vega 


FTS 


1 


5.002 


19.14 


0.02 


AB 


FTS 


V 


9.970 


19.47 


0.03 


Vega 


FTS 


1 


5.475 


19.12 


0.01 


AB 


GS 


V 


16.995 


19.91 


0.04 


Vega 


FTS 


1 


8.972 


19.04 


0.01 


AB 


FTS 


V 


17.899 


19.96 


0.02 


Vega 


FTS 


1 


16.028 


19.12 


0.03 


AB 


FTS 


V 


19.009 


20.03 


0.07 


Vega 


FTS 


1 


17.005 


19.18 


0.02 


AB 


FTS 


V 


23.031 


20.29 


0.04 


Vega 


FTS 


' 


18.958 


19.27 


0.05 


AB 


FTS 


V 


25.939 


20.61 


0.04 


Vega 


FTS 


1 


19.940 


19.34 


0.02 


AB 


FTS 


Y 


35 423 


21 16 


06 


Vega 


FTS 




23 013 


19 55 


03 


AB 


FTS 


K 


0.466 


20.21 


0.01 


Vega 


GS 




25.977 


19.59 


0.05 


AB 


FTS 


K 


1.467 


20.09 


0.01 


Vega 


GS 




33.892 


20.26 


0.08 


AB 


FTS 


R< 


2.028 


20.03 


0.03 


Vega 


FTS 




37.014 


20.63 


0.08 


AB 


FTS 


R c 


2.462 


19.96 


0.01 


Vega 


GS 




46.456 


20.85 


0.07 


AB 


FTS 


R c 


2.891 


19.87 


0.02 


Vega 


FTS 




49.393 


20.99 


0.06 


AB 


FTS 


K 


3.493 


19.77 


0.01 


Vega 


GS 




56.887 


21.20 


0.09 


AB 


FTS 


R< 


5.003 


19.47 


0.02 


Vega 


FTS 




0.470 


19.91 


0.11 


AB 


GS 


K 


5.471 


19.46 


0.01 


Vega 


GS 




3.490 


19.40 


0.01 


AB 


GS 


R c 


7.891 


19.33 


0.02 


Vega 


FTS 




5.470 


19.15 


0.01 


AB 


GS 


R c 


9.072 


19.29 


0.02 


Vega 


FTS 















"Telescope key: FTS: 2m Faulkes Telescope South; GS: 8.1m Gemini-South Telescope. All magnitudes have been corrected for foreground extinction. 



TABLE 3 

Gemini-S Host Photometry of XRF 1003 16D / SN 201 Obh 



Filter 


Mag 


a (mag) 


Telescope 


g 


17.40 


0.08 


Gemini-S 


r 


17.14 


0.07 


Gemini-S 


i 


17.07 


0.06 


Gemini-S 


z 


17.03 


0.06 


Gemini-S 



All magnitudes have been corrected for foreground extinction. 



curves of GRB-SNe 1998bw and 2006aj, as well as local 
SN 19941, all of which have been corrected for foreground 
and host extinction using the values quoted in the literature: 
SN 19941: E(B - V) total = 0.45 mag (Richmond et al. 1996; 
Iwamoto et al. 1996; though see Sauer et al. 2006 who find 
E(B — V) = 0.30 mag from spectral modelling); SN 1998bw: 
E(B - y) tot ai = 0.07 mag (Patat et al. 2001); SN 2006aj: 
E(B — V) to tai = 0.14 mag (Sollerman et al. 2006; though see 
Campana et al. 2006 who find E(B - V) 

host — 0.20 mag from 
their UVOT obsverations.). No other changes have been made 
to the colour curves apart from those of SN 1998bw, which 
have been stretched by s = 0.63 in B — V and V — R c and by 
s = 0.61 in R c — i (see Table q). The stretch factors, which 
are the average values of s determined in the individual filters 
(i.e. s = 0.62,0.64 for respective filters B and V; thus the av- 
erage value is s = 0.63), have been applied so as to present a 
consistent analysis of the properties between SN 201 Obh and 



SN 1998bw (i.e. as the stretch factors were also applied to the 
LCs of SN 1998bw in Figure 

Starting before maximum light, all of the colours of 
SN 201 Obh become redder over time, with similar behaviour 
seen in the comparison SNe. Additionally, SN 201 Obh ap- 
pears to be redder than SN 1998bw and SN 19941 in the 
pre-maximum phases, while in the post-maximum phase 
SN 201 Obh appears to have colours similar to those seen in 
SN 19941. 

3.3. Amis 

We have determined the Amis parameter (i.e. the amount 
that the light curves fades by from peak to fifteen days later) 
of the BVR c i light curves of XRF 100316D/ SN 2010bh. We 
have used low-order polynomials to determine the peak times 
of each light curve, and in turn the Amis parameter. It is seen 
that the light curves decay more slowly in the redder filters, 
with the Amis parameter in filters BVRJ, respectively, being 
1 .80 ± 0. 14, 0.95 ± 0. 10, 0.90 ± 0.05 and 0.77 ± 0.20 (which 
are also listed in Table [5}. The uncertainties in the Ami 5 pa- 
rameters are derived from the uncertainties in determining the 
peak times from the photometry, which, for example, in the 
i band are less certain due to the sparseness of photometry 
around and just after the peak, as well as from the distribu- 
tion of the peak times determined from fitting different order 
polynomials. 

The V- and/? -band values of Am 15 measured for SN 2010bh 
are somewhat larger (i.e. fade faster) than the average values 
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TABLE 4 

Main Photometric parameters ofXRF 1003 16D / SN 2010bh 





B 


V 


Rc 


i 


peak apparent magnitude 


20. 13 ±0.08 


19.47 ±0.08 


19.29 ±0.08 


19.05 ±0.08 


peak absolute magnitude" 


-18.27 ±0.08 


-18.62 ±0.08 


-18.60 ±0.08 


-18.62 ±0.08 


Tpeak (days) 


8.02±0.12 


8.86 ±0.07 


11.18±0.05 


1 1.21 ±0.22 


^v.fore+host (mag) 


1.32 ±0.08 


1.01 ±0.08 


0.81 ±0.08 


0.59 ±0.08 


Am 15 (mag) 


1.80 ±0.14 


0.95 ±0.10 


0.90 ±0.05 


0.77 ±0.20 



" Using a distance modulus of fi = 37.08 mag (for H a = 71km/s/Mpc, Q. M = 0.27 and Q A = 1 - Q. M = 0.73). 
Apparent magnitudes are not corrected for foreground or host extinction. Absolute magnitudes are extinction corrected (foreground & host). 
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FIG. 3.— B-band light curve of XRF 100316D / SN 2010bh (solid points), 
which has been corrected for foreground and host extinction. Plotted for com- 
parison are the B-band light curves of two known GRB-SNe: 1998bw and 
2006aj. The light curves of the comparison SNe have been shifted by dif- 
ferent amounts in each filter to match in peak brightness. No other changes 
have been made to the light curves apart for the light curves of SN 1998bw 
which have been stretched by a factor s = 0.62 (see Table |5J. It is seen 
that SN 2010bh did not increase linearly in brightness from first detection 
like SN 1998bw, and at t — 1„ = 0.59 days is of comparable brightness as the 
shifted light curve of SN 2006aj. 

measured for a sample of local Ibc SNe observed in the V 
band and R band by Drout et al. (201 1) (their Table 4). It is 
also seen that SN 2010bh fades faster, on average, than la SNe 
in the B-band (e.g. Phillips 1993; Riess et al. 1998). 

3.4. Shock Break-out or optical afterglow? 

The break-out of the blast/shock-wave and/or emission 
coming from the shock-heated stellar envelope, has been de- 
tected in both type Ibc and type II SNe, including: type Hp 
SN 1987A (Early peaks in optical LCs: Hamuy et al. 1988, 
and PCyg-like feature around 1500 A in UV light echo off 
of a dust cloud ~ 300 pc from the SN: Gilmozzi & Panagia 
1999); type lib SN 1993J (Early peak in optical light curves: 
Wheeler et al. 1993, Schmidt et al. 1993, Richmond et al. 
1994, Lewis et al. 1994); type Ibc SN 1999ex (early optical 
peaks: Stritzinger et al. 2002); type IIP SNe SNLS-04D2dc 
& SNLS-06Dljd (UV-flash: Gezari et al. 2008); type IIP 
GALEX supernova SNLS-04D2dc (UV-flash: Schawinski et 
al. 2008); type Iln PTF 09UJ (UV-flash: Ofek et al. 2010); 
and type IIP SN 2010aq (early UV & optical peaks: Gezari et 
al. 2010). 

One SN in particular has been the subject of much inter- 
est and debate when intrepreting the high energy properties 
of the event. SN 2008D (Modjaz et al. 2008; Soderberg et 
al. 2008; Mazzali et al. 2008; Malesani et al. 2009) oc- 




0.01 



10000 
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FIG. 4.— Spectral Energy Distributions (SEDs) of XRF 1003 16D / 
SN 2010bh over the first 8.6 days. All photometry has been corrected for 
foreground and host extinction. The first five epochs (t — 1„ ks 0.5, 1.5, 
2.5, 3.5 and 5.5 days) are from the Gemini-S detections, while the SED 
at t — 1„ = 8.5 days is from the first HST data. The power-law index (i.e. 
F v oc v?) of the first SED is j3 = +0.94 ±0.05; £ 2 / do - f = 2.17/2 (best- 
fitting line not plotted). The slope of the spectrum is steeper than that ex- 
pected for light emitted as synchrotron radiation, where, when neglecting 
effects due to self-absorption, a maximum value of j3 = +1/3 is allowed in 
synchrotron theory (e.g. Sari et al. 1998). Thus, the origin of light is not ex- 
pected to be synchrotron but rather is coming from the shock break-out from 
the stellar surface and dense stellar wind. 



curred in nearby galaxy NGC 2770 while Swift was observ- 
ing SN 2007uy within the same galaxy. The resultant high- 
energy emission detected by Swift was interpreted differently 
in the literature, with the hot, black-body X-ray spectrum & 
early peak in optical LC being intrepreted by Soderberg et al. 
2008 as being due to the shock break-out, while Mazzali et 
al. (2008) suggests that the emission is due to a "choked" jet. 
Interestingly, a recent paper by Van der Horst et al. (2011) 
has perhaps put the debate to rest. Using observations of the 
radio emission of SN 2008D during the first year after the 
explosion, the authors showed that there was no evidence of 
a relativistic jet contributing to the observed radio flux, thus 
strongly suggesting that the high-energy emission was due to 
the shock break-out and not to a GRB- or XRF-like event. 

Emission coming from the shock-heat, expanding stellar 
evelope as also been detected in GRB/XRF-SNe, though the 
topic is still contentiously debated in the literature. For XRF 
060218, the hot, black-body X-ray spectrum & early peak in 
UV and optical LCs was intrepreted by Campana et al. 2006 
(C06 here-on), and later Waxman et al. (2007), as being due to 
the shock break-out. However, the emission was intrepreted 
by Ghisellini et al. (2006) as being due to a relativistic jet. 

Figure [3] shows the B-band LC of SN 2010bh, which has 
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been corrected for foreground and host extinction. Plotted for 
comparison are the light curves of two GRB-SNe: 1998bw 
and 2006aj, which have been shifted to match SN 2010bh 
in peak brightness, and the LC of SN 1998bw has also been 
stretched by a factor s = 0.62 (see Table[5]l. 

Upon inspection, similarities between SN 2010bh and 
SN 2006aj are seen in the early-time LC, where a bright 
component is detected at t — 1 = 0.598 days. The magni- 
tude at this epoch is similar to that at t — t Q = 3.489 days, in- 
dicating that SN 2010bh did not increase linearly in bright- 
ness from detection, as was seen for SN 1998bw. In the 
case of SN 2006aj, early peaks in the UV and optical LCs 
(t — 1 < 10 4 s) were explained by C06 as light coming from 
the low energy tail of the thermal X-ray emission (see below) 
produced by the radiation shock driven into the stellar wind. 
At later times, the optical and UV emission is attributed to the 
expanding envelope of the progenitor star that is heated as the 
shock-wave passes through it. Initially the envelope is opaque 
due to the dense stellar wind, but as the star and wind expands, 
the photosphere propogates inward exposing the shocked stel- 
lar material. For SN 2010bh, though we are not able to ascer- 
tain the exact shape of the early Z?-band light curve due to the 
paucity of our early-time observations, it is not unreasonable 
to envision a scenario similar to that observed for SN 2006aj 
in which the brightness of the flux at this early epoch is due 
to flux coming from the shock-heated, expanding stellar en- 
velope. 

In addition to the bright B-band detection, inspection of the 
spectral energy distribution (SED) at this epoch, which has 
been corrected for foreground and host extinction, gives ad- 
ditional clues to the origin of the flux. At the first epoch the 
SED is very blue and does not resemble the shape of the SEDs 
at later times (see Fig|4|i. 

To determine whether the flux at this time is synchrotron 
in origin, the optical SED was fit with a power-law (i.e. 
F v oc vP), for which we find j8 = +0.94 ± 0.05 (* 2 /d.o.f = 
2.17/2). The slope of the spectrum is harder than that ex- 
pected for light emitted as synchrotron radiation, where, when 
neglecting effects due to self-absorption (see below), a max- 
imum value of /3 = +1/3 is allowed (e.g. Sari et al. 1998). 
Thus, the shape of the SED at this early epoch is consistent 
with emission coming from the shock-heated stellar envelope. 

Additionally, we have constructed SEDs for SN 2006aj at 
t — 1 = 0.50 days and 0.75 days and fit them with a sin- 
gle power-law. Quite interestingly, at t — t Q = 0.50 days we 
find j8 = +0.97 ± 0.05 (% 2 /d.o.f = 2.13/2), while at t - 1 = 
0.75 days we find j3 = +0.94 + 0.03 (x 2 /d.o.f = 2.03/2). 
These values for the spectral index are akin to that found for 
SN 2010bh at a similar epoch, and are both harder than that 
expected for synchrotron radiation. 

The assumption that the optical flux is above the self- 
absorption frequency (v„) cannot be proved with our current 
set of observations, however it has been seen in other GRB 
events that the synchrotron self-absorption frequency at typi- 
cal observing epochs of about 1 day after the GRB is expected 
to be of order ~ 10 9 — 10 10 Hz (i.e. at radio frequencies). In 
the literature there are several events where there has been 
enough multi-wavelength observations to determine the self- 
absorption frequency, including: GRB 970508: V„~3x 10 9 
Hz (Granot et al. 1999); V a - 2.5 x 10 9 Hz (Galama et al. 
1998); GRB 980329: v„ - 13 x 10 9 Hz (Taylor et al. 1998); 
GRB 991208 v„ ~ 4 - 1 1 x 10 9 Hz (Galama et al. 2000). As 
the frequency range of the optical observations lie between 



^ 3 — 10 x 10 14 Hz (i.e. ~ 10 4 — 10 5 times higher frequency), 
it is unlikely that the optical observations suffer from self- 
absorption effects. 

Further evidence that the light at t — 1 = 0.598 days is 
coming from the shock-heated stellar envelope is the analy- 
sis performed by Starling et al. (2011) on the X-ray spec- 
trum at 144 s < t — 1 < 737 s. The authors found evidence 
for a soft, hot, black-body component that contributes « 3% 
to the total 0.3 — 10.0 keV flux (though see also Fan et al. 
2011 for a different interpretation). With this in mind, the 
presence of the black-body component, with kT = 0.14 keV 
(T r; 1.6 x 10 6 K) is similar in temperature to the soft, hot, 
black-body component seen in the X-ray spectrum of XRF 
060218 / SN 2006aj (kT w 0.17 keV; T ~ 2.0 x 10 6 K; C06). 

According to C06, the thermal component is key to under- 
standing XRF 060218. It is known that the signature of shock 
break-out is a hot black-body X-ray spectrum immediately af- 
ter the explosion, and C06 argued that the hot temperature of 
the thermal X-ray component is indicative of radiation emit- 
ted by a shock-heated plasma. Furthermore, the characteristic 
radius of the emitting region was ~5x 10 12 cm (Waxman 
el al. 2007 also find the emitting region to be at a radius of 
7.8 x 10 12 cm), which though is quite large, (i.e. of order of 
the radius of a blue supergiant), was explained by C06 and 
Waxman et al. (2007) by the presence of a massive, dense 
stellar wind surrounding the progenitor, which is common for 
Wolf-Rayet stars (the supposed progenitors of type Ibc SNe 
and possibly long-duration GRBs). Both authors argue that 
the thermal radiation is observed once the shock, that is driven 
into the wind, reaches a radius where the wind becomes opti- 
cally thin. 

For XRF 100316D / SN 2010bh, the temperature of the 
black-body component in the X-ray spectrum implies an emit- 
ting radius of ~ 8 x 10 11 cm, which is almost an order of 
magnitude larger than the radius of a Wolf-Rayet star (~ 10 11 
cm). As for XRF 060218, the large radius of the X-ray emit- 
ting region could be the result of the shock being driven into 
the dense stellar wind, and the thermal radiation is observed 
once the wind density decreases and becomes optically thin. 

Thus, the brightness of the B-band detection at t — t Q = 
0.598 days, when taken in tandem with the hard value of 
the power-law index (J3 = +0.94 + 0.05) of the SED at this 
epoch, as well as the presence of the extremely hot thermal 
component (T ps 1.6 x 10 6 K) in the X-ray spectrum at 144 s 
< t — to < 737 s, suggests that we have detected optical emis- 
sion coming from the cooling, expanding envelope that was 
heated by the passage of the shock-wave through it. 

3.5. Stretch Factor relative to SN 1998bw 

We have determined the stretch factor (s) and luminosity 
factor (k) of SN 2010bh in relation to the archetype GRB- 
SN, SN 1998bw. To do this we have created synthetic flux 
SN 1998bw light curves in filters BVR c i as they would ap- 
pear if they occurred at z = 0.0591. At each epoch, an SED 
of the observed SN 1998bw light curves is built and interpo- 
lated using polynomials. This allows us to cover the temporal- 
frequency space for SN 1998bw and build synthetic light cu- 
vres for other specified observed frequencies. When con- 
structing the synthetic LCs we have accounted for the inverse- 
squared distance (i.e. luminosity distance) suppression of 
flux, as well as rest-frame extinction and time dilation effects. 
We then fit the flux SN 1998bw light curves with an empirical 
relation of the form: 
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TABLE 5 

Luminosity (k) and Stretch (s) Factors of SN 201()bh 



Filter 


s 


k 


B 


0.62 ±0.01 


0.41 ±0.01 


V 


0.64 ±0.01 


0.43 ±0.01 




0.62 ±0.02 


0.40 ±0.01 


i 


0.60 ±0.01 


0.48±0.01 



U(t)=A + st[ e — - +r g Mog(r g 3) (1) 

where A is the intercept of the line and s and F are related 
to the respective amplitude and width of the function. The 
exponential cut-off function for the rise has a characteristic 
time R and a phase zero-point p, while a\ , OC2 and 0C3 are free 
parameters. All of the parameters are allowed to vary during 
the fit. 

Once the flux light curve of SN 1998bw was fit by Equation 
(Q}, the stretch and luminosity factors of SN 2010bh were de- 
termined by fitting the following equation: 



W(t)=kxU{t/s) (2) 

to the flux light curve of SN 2010bh in each filter. 

The stretch factors of SN 2010bh relative to SN 1998bw 
were determined using data taken several days after the ini- 
tial burst (t — 1 > 7.0 days) as the early light curves, whose 
shapes are altered by light coming from an additional shock 
break-out component, are very different to that of SN 1998bw. 
The results of the fit are listed in Table [5] where the quoted er- 
rors are statistical only. 

For GRB and XRF events where an optically bright SN has 
been detected (i.e. not including events such as SN-less GRBs 
060505 and 060614; e.g. Fynbo et al. 2006), it is seen that 
SN 2010bh is the faintest SN to date that has been linked 
spectroscopically with a GRB or XRF (see Table |HJ, and is 
of similar peak brightness as the SN associated with GRB 
970228 (k = 0.40 ± 0.29, though the host extinction is un- 
known), but not as faint as the (possible) supernova that was 
photometrically-linked to GRB 101225A (k = 0.08 ±0.03, 
Thone et al. 2011). The slower evolution of SN 2006aj rel- 
ative to SN 2010bh seen upon inspection of the optical light 
curves are echoed in the computed stretch factors for the two 
SNe. In the B filters they are approximately the same value 
(SN2010bh: s — 0.62±0.01; SN 2006aj: s = 0.60 ±0.01), 
however in the redder filters SN 2010bh has smaller stretch 
factors relative to SN 2006aj. 

It should be mentioned that SN 1998bw is not the ideal 
template for SN 2010bh. At early times light from the shock 
break-out changes the shape of the light curve, which affects 
how well SN 1998bw, for which no shock break-out was ob- 
served, can be used as a template. Even at late times it is 
seen in Figure[T]that SN 2010bh decays more slowly than the 
stretched LC of SN 1998bw, thus showing that the two SNe 
evolve quite differently. The limitations of using SN 1998bw 
as template was also noted by Ferrero et al. (2006), where 
they used an additional power-law component when fitting 
their SN 1998bw template to their multi-filter observations of 
SN 2006aj. 



4. REDDENING & ABSOLUTE MAGNITUDES 

Determining the absolute magnitude of SN 20 lObh requires 
an estimate of the amount of extinction local to the site of the 
SN explosion. Drout et al. (2010) found for a sample of local 
type Ibc SNe that the mean values of the (V — R) colour at ten 
days after maximum are tightly distributed. The authors found 
at ten days after V maximum ((V —R)vw) = 0.26 ±0.06 mag, 
and ten days after R maximum ((V -R)mo) = 0.29 ±0.08 
mag. For SN 2010bh we found the (V —R) colour at ten days 
after V and R maximum to be, respectively, ((V — R)viq) = 
0.38 ± 0.05 mag, and ((V - R) RW ) = 0.43 ± 0.05 mag. The 
difference in the colours of the sample found by Drout et al. 

(2010) in comparison with SN 201bh implies a colour excess 
of E(V — R) ~ 0.12 mag, which in turn implies E(B — V) ~ 
0.18 mag (using the relative extinction values in Table 6 of 
Schlegel et al. 1998). 

For comparison, using the Balmer Ha and Hj3 lines fluxes 
to measure the Balmer decrement, Starling et al (201 1) found 
a combined reddening (host and foreground) of E(B — V) = 
0.178 mag for a nearby HII region ("Source A"; see Fig- 
ure 11 in Starling et al. 2011). The foreground extinc- 
tion found from the dust maps of Schlegel et al (1998) is 
E(B — V) = 0.1 17 mag, which implies a rest-frame, host ex- 
tinction of E(B — V) = 0.061 mag. This value is consistent 
with the value found via comparing the colour curves, al- 
beit somewhat smaller, but provides additional credence to 
the method described by Drout et al. (2010). We note that 
the value for the colour excess found by Starling et al. (201 1) 
was not for the specific explosion site of SN 2010bh as the 
SN was still bright during their observations, thus the au- 
thors made the colour excess measurement of a nearby bright 
knot/blob instead. None the less, these independent analysis 
imply a small colour excess for the location near and around 
SN 2010bh. Using the value E(B-V) =0.18 ±0.08 mag, 
we have calculated the absolute magnitudes for SN 2010bh 
(Table!). 

In comparison with the peak magnitudes of almost all of 
the previously detected GRB-SN tabulated by Cano et al. 

(2011) , it is seen that the peak absolute V-band magnitude 
of SN 2010bh {My = -18.62 ±0.08) is the faintest ever de- 
tected for a spectroscopically-detected GRB-SNe and is ~ 0.8 
mag fainter than SN 1998bw at peak brightness. Cano et al. 
(2011) found for a sample of 22 GRB-SNe (the authors ex- 
cluded events where only an upper limit to the SN brightness 
was obtained; e.g. GRBs 060505 & 060614, Fynbo et al. 
2006), where measurements of the host extinction have been 
made and included in the absolute magnitude calculation, an 
average, absolute, V-band magnitude of (My) = —19.00, with 
a standard deviation of a = 0.77 mag. In a separate analysis 
of the peak GRB-SNe magnitudes, Richardson (2009) deter- 
mined for a sample of 14 GRB-SNe with known values for 
the host extinction, an average absolute V-band magnitude of 
(My) = —19.20, with a standard deviation of a = 0.70 mag. 
In comparison with these analyses it appears that SN 2010bh 
lies at the lower end of the peak brightness distribution. 

5. EXPLOSION PARAMETERS 

We have constructed a quasi-bolometric light curve of XRF 
100316D / SN 2010bh (Figure EJ in the 3, 000 A - 16, 600 A 
wavelength range (i.e. UBVR C UH) by taking the follow- 
ing steps: (1) converting all magnitudes to monochromatic 
fluxes after correcting for foreground and host extinction; (2) 
the SED at each epoch was integrated over frequency. At 
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Time Since Burst (days) 

FIG. 5.— The quasi-bolometric light curve for XRF 100316D/SN 2010bh 
in the 3,000A— 16, 60()A wavelength range, which was obtained by integrat- 
ing the flux in the optical and infrared filters UBVRiJH. The model (solid 
line) that is used to determine the nickel mass and ejecta mass is also plotted. 
Peak bolometric light is found to be at ; — 1 = 8.57 ± 0.04 days. 



applicability of the diffusion approximation for photons (i.e. 
the photospheric phase). We also assume a constant opacity 
K = 0.07 cm 2 g _1 (e.g. Chugai 2000), which is justified if 
electron scattering is the dominant opacity source0 (e.g. 
Chevalier 1992). Finally, with respect to Arnett (1982) 
who only considered the energy produced by the decay of 
nickel into cobalt ( 56 Ni — > 56 Co), we also include, using the 
methodology described in Valenti et al. (2008), the energy 
produced by the decay of cobalt into iron ( 56 Co -> 56 Fe): 

L(t) = M Ni e~ x2 x 



(^(eN 1 -eco)^ X A(z)dz + eco^ X B(z)dzj (3) 

where 



epochs where there was no contemporaneous data, the indi- 
vidual light curves in each filter were extrapolated after fitting 
Equation (fl~|i to the individual light curves in each filter. For 
epochs where there was no IR data, we estimated the amount 
of IR flux as a fraction of the total bolometric flux using the 
existing data and extrapolating between epochs. 

For epochs before the first IR detection (i.e. t — t c = 8.6 
days), we have assumed a constant fraction of flux at IR wave- 
lengths of 20%. We note that this assumption of a constant 
amount of flux at IR wavelengths is valid only for the first 
dozen or so days, as it is seen that the fraction of light at IR 
wavelengths actually increases with time and can be as high 
as 45% at late times (see Figure|8]in this work, as well as Fig- 
ure 9 of Modjaz et al. 2009). The increase in the amount of 
flux at redder wavelengths is expected however, as it is seen 
that core-collapse SNe become redder over time. We estimate 
the uncertainty in the bolometric magnitude, which is domi- 
nated by systematic errors described in the method above, to 
be of order 0.2 mag. 

We have calculated the time of peak bolometric light by 
fitting a series of different order polynomials to the bolo- 
metric light curve. We find peak bolometric light to be at 
t — 1 = 8.57 ±0.04 days, where the error is statistical and de- 
termined from the distribution of peak times determined from 
the different order polynomials. 

Here we have determined a quasi-bolometric light curve of 
XRF 1003 16D / SN 2010bh up to re 35 days. This time inter- 
val corresponds to the photospheric phase. To attempt to ex- 
tract physical parameters of the progenitor star we have used a 
simple analytical model that was developed by Arnett (1982) 
for type I SNe, which was later expanded by Valenti et al. 
(2008). 

Arnett-like models have been used by many authors to 
intrepret the light curves of numerous Ibc SNe and GRB-SNe 
events, including: lb SN 1999dn, Benetti et al. (2011); Ic-BL 
SN 2003jd, Valenti et al. (2008); Ic SN 2004aw, Tauben- 
berger et al. (2006); Ic-BL SN 2009bb, Pignata et al. (201 1); 
GRB-SNe 1998bw, Iwamoto et al. (1998); Richardson et 
al. (2006) for a sample of Ibc SNe; Drout et al. (2011) for 
a sample of Ibc SNe. The model assumes: (1) homologous 
expansion of the ejecta, (2) spherical symmetry, (3) all of the 
radioactive nickel ( 56 Ni) is located at the centre of the explo- 
sion and does not mix, (4) radiation-pressure dominance, (5) 
a small initial radius before explosion (R — > 0), and (6) the 



A(z) =2ze- 2zy+z2 ,B{z) = 2ze 



-2zy+2zs+z 2 



(4) 



and x ee t/x m , y = x m /(2x m ), and s = (x m (x Co ~ 
TiV()/(2Tc„TA?,)). The energy release in one second by one 
gram of 56 Ni and 56 Co are, respectively, e Ni = 3.90 x 10 10 
erg s _1 g _1 and Eco = 6.78 x 10 9 erg s _1 g _1 (Sutherland 
& Wheeler 1984; Cappellaro et al. 1997). The decay times 
of 56 Ni and 56 Co, respectively, are Xm — 8.77 days (Tauben- 
beger et al. 2006 and references therein) and Xq = 1 1 1 - 3 days 
(Martin 1987). 

X m is the effective diffusion time and determines the width 
of the bolometric light curve. x m is expressed in relation to 
the opacity k and the ejecta mass M e j, as well as the photo- 
spheric velocity v p h at the time of bolometric maximum and 
the kinetic energy of the ejecta E^: 



K 



V2/M2A 1/4 



(5) 



where /3 re 13.8 is a constant of integration (Arnett 1982). 
Thus, for SNe with similar ejecta velocities, a wider bolomet- 
ric light curve implies more ejecta mass and higher kinetic 
energies, while for SNe with similar amount of ejected mate- 
rial, a wider light curve implies slower ejecta velocities. 

We have fit the data at times > 5 days (as to not include 
flux that may also be coming from the shock-heated, ex- 
panding stellar envelope). We have taken the photospheric 
velocity around peak light for XRF 100316D / SN 2010bh 
to be Vph re 25,000 km s (determined from the Si II ab- 
sorption lines measured by Chornock et al. 2010; their Fig- 
ure 3), the analytical model yields M Ni = 0.10 ± 0.01M K „ 
M ej = 2.24 ± 0.08M.. E k = 1 .39 ± 0.06 x 10 52 erg. Our re- 
sults are displayed in Table [6] where they are compared with 
other Ibc and GRB-SNe. 

The quoted errors are statistical and include the uncertain- 
ties in the photometry, extinction, and calibration. However, 

22 Metal lines also contribute to the total opacity, and strictly speaking the 
electron scattering opacity is not actually constant in time as the ionization 
state changes. For triply ionized ejecta, K R3 0.07. 
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the dominant sources of error are systematic. For example, 
our assumption of a constant value of the opacity, is perhaps 
valid during the photospheric phase, but may differ in value 
from our assumed value of K = 0.07 cm 2 g _1 . Using Equa- 
tion [5] and the values of the explosion energy and ejecta mass 
from the literature for SN 1998bw, we calculate an opacity 
of K « 0.05 cm 2 g _1 . Since the total opacity depends on the 
metal content of the ejecta as well as that from electron scat- 
tering, it is quite possible that the value of the opacity may 
have a slightly higher or lower value than what we have used 
here. A small change in the value of the opacity leads to an 
uncertainity in the ejected mass of (M e j) of 25% — 30% and 
an uncertainity in the kinetic energy (E^) of 60% — 75%. The 
uncertainty in the nickel mass is less, as it is mainly affected 
by the uncertainty in the rest-frame extinction, and to a lesser 
extent in the photometric errors. 

As an alternative method to check how much mass was 
ejected during SN 2010bh, we have compared the width 
of bolometric LC of XRF 100316D / SN 2010bh with 
those of three other spectroscopically-confirmed GRB-SNe: 
SN 1998bw, SN 2003dh and SN 2006aj. If we assume that 
the opacity is the same in each event, we can solve Equation 
(01 to find the ejected mass for SN 2010bh: 
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FIG. 6.— The UVOIR bolometric light curves, in the 1,900A- 16,660A 
wavelength range, of XRF 100316D / SN 2010bh (filled circles) and XRF 
060218 / SN 2006aj (grey hexagons). When constructing the bolometric 
light curve of SN 2010bh, we have assumed that the same fraction of the 
bolometric flux emitted at UV wavelengths is the same for SN 2010bh and 
for SN 2006aj. Plotted for comparison is the UBVR C UH bolometric light 
curve of SN 2010bh (open circles). All photometry has been corrected for 
foreground and host extinction. 



/vi0bh\ /Tm,10bh\ 

Mej.lObh «M ej .GRB (6) 

\VgRB/ \T m ,GRB/ 

First, using SN 1998bw, from the literature (e.g. Iwamoto 
et al. 1998; Patat et al. 2001) it is expected that m 8 - 10M@ 
of material was ejected. We have measured Z m psbw fr° m our 
bolometric LC of SN 1998bw (see Section 7.2, where we 
have also measured % m for SN 2006aj and SN 2003dh on way 
to deriving physical parameters using our analytical model) 
and thus derive M e j 10 bh ~ 3.4M . For SN 2003dh, it is ex- 
pected that « 7 ± 3M S of material was ejected (e.g. Deng et 
al. 2005), which in turn implies for SN 2010bh an ejected 
mass of Mej.iobh ~ 3.3M©. For SN 2006aj, it was deter- 
mined that w 2M© of material was ejected during the explo- 
sion (e.g. Mazzali et al. 2006). This implies an ejected mass 
for SN 2010bh of M e j. 10 bh ~ 2.7M Q . In turn, these values for 
the ejected mass imply, respectively, explosion energies for 
SN 2010bh of E kil0bh w 2.1 x 10 52 erg, E k . 10bh « 2.0 x 10 52 
erg, and E k iobh ~ 1.7 X 10 52 erg. This method gives larger 
values for the ejected mass of SN 2010bh which imply larger 
explosion energies. On average, these values are larger than 
those determined directly from our fit. 

It is important to consider the limitations of our model 
and the implied caveats from our assumptions. While it 
appears that our model is able to accurately reproduce the 
nickel masses in each of these events, the energy estimates 
have probably the largest uncertainities. As mentioned, Ar- 
nett (1982) considers spherical geometry, while GRB-SNe are 
thought to be asymmetric and may well have an unusual dis- 
tribution of nickel mass in the ejecta mass. Furthermore, this 
model calculates an explosion energy by assuming that all of 
the ejecta are moving at the photospheric velocity. In the case 
of XRF 1003 16D / SN 2010bh, while Chornock et al. (2010) 
observed absorptions in the spectra at 25,000 km s 1 , it does 
not mean that all of the ejecta are moving at that velocity. It 
is quite possible that only a small or modest fraction of the 
mass is expelled at high velocities (either isotropically or col- 
limated into a jet) creating the absorption lines in the spectra, 
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FIG. 7.— XRF 060218 / SN 2006aj: The amount of UV flux as a fraction 
of the total bolometric flux in the 1,900A — 16,660A wavelength regime, 
plotted over time. It is seen that up to a few days after the initial explosion, 
the amount of UV flux is considerable. More importantly, even at late times 
(> 10 days) the amount of flux at UV wavelengths is a non-negligible fraction 
of the total bolometric flux. 

while most of the mass moves at slower velocities. A more ac- 
curate way of determining the actual explosion energy would 
be to understand the way the mass is distributed in the ejecta 
(i.e. via detailed modeling), and integrate over the density- 
velocity profile. 

6. DISCUSSION AND FURTHER ANALYSIS 
6.1. SN 201 Obh in relation to SN 2006aj 

It was shown in the preceeding section that, despite ejecting 
only a moderate amount of material, SN 201 Obh was an ener- 
getic explosion, mainly due to the large photospheric veloci- 
ties. To put the results of our bolometric light curve modelling 
into context, we first compare SN 201 Obh with SN 2006aj and 
then with the other GRB-SNe. 

The ultra-violet, optical and infrared (UVOIR) bolometric 
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FIG. 8. — The amount of flux as a fraction of the total bolometric flux (in 
filters UBVR C UH) that is obtained by integrating the optical (grey, dashed- 
line) and infrared (black, solid-line) flux for SN 2()10bh and SN 2009bb. It is 
seen that the behaviour is similar in each SN event. 



light curve (in the 1,900 A — 16, 660 A wavelength regime) of 
XRF 060218 / SN 2006aj, which has been constructed from 
the data published in the literature (Campana et al. 2006; 
Sollerman et al. 2006; Cobb et al. 2006; Ferrero et al. 2006), 
is displayed in Figure|6] It is seen that during the first few days 
the fractional contribution of UV flux (i.e. light detected in the 
1 , 900 A — 3, 000 A wavelength regime) to the total bolometric 
LC (Figure |7]i is considerable due to the shock break-out of 
the progenitor star (Campana et al. 2006), which decays at a 
linear rate after « 4 days. It is important to note that the con- 
tribution of flux at UV wavelengths is not negligible, and at 
roughly 10 days still contributes 15% of the total bolometric 
flux. 

By making the assumption that the fraction of UV flux 
of the total bolometric flux at a given epoch is the same in 
SN 2010bh as it is in SN 2006aj, we have used SN 2006aj 
as a template/proxy for creating a "UV-corrected" bolometric 
light curve for SN 2010bh, which is also displayed in Figure 
[6] The bolometric light curve of SN 2010bh is shown twice: 
(1) first, by considering flux only from the optical (UBVRJ) 
and infrared (JH) filters and, (2) by including a contribution 
of flux at UV wavelengths. 

Our assumption that the fraction of UV flux is similar in 
each event is perhaps not unreasonable when we also consider 
the percentage of the bolometric flux detected in the optical 
and IR regimes for SN 2010bh and SN 2009bb. Plotted in 
Figure[8]is a plot of the amount of light detected in the optical 
and infrared as a fraction of the total bolometric flux in each 
epoch for SN 2009bb and SN 2010bh. The data included in 
this plot does not include all of the data obtained for these 
two events, but only those where there is contemporaneous 
optical and IR data. The first IR detection of SN 2010bh was 
during the first HST epoch at t — t G = 8.6 days, and the first 
IR detection of SN 2009bb was at t - to = 7.7 days (Pignata 
et al. 201 1). Inspection of the optical and infrared filters pairs 
reveal similar behaviour in both events, though it is seen that 
more flux is emitted at IR wavelengths in SN 2010bh relative 
to SN 2009bb. 

The fractional contribution of the optical and IR flux to the 
bolometric light curves for these two events, as well as the 
large UV contribution seen in SN 2006aj, is similar to that 
seen in XRF 080109 / SN 2008D. Modjaz et al. (2009) also 



compared the relative contribution of the UV, optical, and IR 
flux in the bolometric LC of SN 2008D for the first w 30 days 
(Figure 9; Modjaz et al. 2009), and noted a large contribution 
of flux from 2008D at UV wavelengths (» 40% at t- to « 0.8 
days). This is not quite as high as for SN 2006aj which is 
w 70% at t — 1 w 0.8 days, but their UV contribution does 
not include far-UV flux, whereas SN 2006aj does. Indeed 
at t — 1 = 1.8 days, when they do include the far-UV flux, 
the total UV flux is seen to increase. Modjaz et al. (2009) 
also showed that the amount of IR flux increased over time, 
while the optical flux increased linearly for the first ten days 
or so, then decreased. While there is scant data plotted in 
Figure [8] before t — 1 = 10 days, our plot does show that the 
optical contribution decreases during 8.6 days < t — t G < 24.8 
days for SN 2010bh, whereas for SN 2009bb, decreases until 
~ 25 — 30 days, then appears to contribute a nearly constant, 
if slightly increasing, amount to the total bolometric flux. In- 
deed, when comparing Figure [8] with Figure 9 of Modjaz et 
al. (2009), the general behaviour of SNe 2008D, 2009bb and 
2010bh appear to follow similar trends in their relative bolo- 
metric fluxes during the first w 30 days. 

Next, we applied our analytical model to our constructed 
bolometric light curve of XRF 060218 / SN 2006aj. Tak- 
ing the peak photospheric velocity to be v p h rs 20,000 km 
s" 1 (Pian et al. 2006), we find: M Ni = 0.16 ± 0.01M Q , 



Me, = 1 .63 : 0.05M . . E k = 6.47 ±0.20 x 10 M erg. In com- 
parison, Mazzali et al. (2006) found when modeling the spec- 
tra and photometry of SN 2006aj parameters of M^ ~ O.2M , 
M e j « 2M©, Ek w 2 x 10 51 erg. It appears that our simple 
model is consistent with the values of the physical parameters 
of the SN 2006aj explosion determined via detailed modeling, 
but over-estimates the explosion energy by a factor of ~ 3. As 
mentioned in Section (5), our estimate of the energy assumes 
that all of the mass moves at the photospheric velocity. Here is 
a clear example where this assumption over-estimates the ex- 
plosion energy, implying that a bulk of the ejecta are moving 
at slower velocities. 

Alternatively, if we exclude the contribution of flux at UV 
wavelengths, similar to the UBVR C I C JH bolometric LCs pre- 
sented by Pian et al. (2006), we find for SN 2006aj: M Ni = 
O.14±O.O1M , M ej = 1.92±O.12M , E k = 7.64 ± 0.49 x 
10 51 erg, which is again consistent with those found by Maz- 
zali et al. (2006), but the explosion energy is again over- 
estimated, this time by a factor of ~ 4. In both cases the errors 
are estimated as per Section (5). We note that by including the 
contribution of flux at UV wavelengths, which is considerable 
at early times (i.e. < 5 days), this increases the total lumi- 
nosity of the SN, but makes the bolometric light curve peak 
earlier, thus decreasing the overall width of the LC (i.e. less 
mass is ejected). The results of our modelling are displayed 
in Table [6] 

When we apply our model to the bolometric light curve of 
2010bh that includes a UV correction, we find M^i = 0.12 ± 
0.01M Q , M ej = 1.93±0.07M@, E k = 1.20 ±0.05 x 10 52 erg. 
The inclusion of the UV flux makes the overall bolometric 
light curve brighter, implying more 56 Ni was created during 
the explosion. However, as for SN 2006aj, the inclusion of the 
UV flux, which is dominant at early times, causes the light 
curve to peak earlier, thus making the overall shape of the 
bolometric light curve narrower, and implying that less mass 
was ejected during the explosion (and thus a less energetic 
explosion). 

Never the less, whether an estimated contribution of the UV 
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flux is included or not, the results of bolometric LC modeling 
imply that SN 2010bh ejected a mass of M e j = 1.9 - 2.2M , 
has a nickel mass of = 0.10 — 0.12M Q , and a kinetic en- 
ergy of roughly Ek = 1 .2 — 1 .4 x 10 52 erg. 

6.2. Bolometric Light Curves of other GRB-SNe and local 
stripped-envelope SNe 

It is seen that, despite the similaries of the high-energy 
properties of XRF 100316D and XRF 060218 (Starling et 
al. 2011), the accompaning supernova are somewhat differ- 
ent, with SN 2010bh being much more energetic. In this 
section we compare the bolometric properties of SN 2010bh 
with other GRB-SNe, as well as non-GRB-SNe: Ic SN 19941 
and broad-lined Ic SN 2009bb. We investigate whether clear 
differences exist in the bolometric properties of stripped - 
envelope, core-collapse SNe that produce a GRB and those 
that do not. 

To date, numerous authors have determined physical pa- 
rameters of Ibc SNe by modeling photometry and spec- 
troscopy. Many authors employ analytical models that are 
derived from the original prescription of Arnett (1982) to de- 
scribe bolometric LCs, however very often these models only 
supplement the derivation of physical parameters that are de- 
termined by comparing synthetic spectra created by synthesis 
codes to observations. In this work we are attempting to ob- 
tain physical parameters of previously-studied events by mod- 
eling only the bolometric LCs, and compare our results with 
those obtained via more detailed analyses. 

Figure [9] shows the optical & infrared bolometric light 
curves of XRF 100316D / SN 2010bh in relation to three 
spectroscopically-confirmed GRB-SNe: 1998bw, 2003dh & 
2006aj, as well as broad-lined Ic SN 2009bb and Ic SN 19941, 
neither of which were accompanied by a GRB -trigger. All of 
the bolometric light curves have been determined by integrat- 
ing the flux in the optical and infrared UBVR C I C JH filters, and 
have been corrected for foreground and host extinction. We 
have constructed the bolometric light curves of SN 1998bw 
(Galama et al. 1998a; McKenzie & Schaefer 1999; Sollerman 
et al. 2002; Patat et al. 2001), SN 2006aj (Campana et al. 
2006; Sollerman et al. 2006; Cobb et al. 2006; Ferrero et al. 
2006) & SN 2009bb (Pignatta et al. 201 1) from the individual 
light curves presented in the literature. We have also incorpo- 
rated the bolometric data presented in Deng et al. (2005) for 
SN 2003dh and the bolometric data presented in Richmond et 
al. (1996) for SN 19941. When constructing the bolometric 
light curve of SN 1998bw, for which very little infrared data 
are available in the literature, we estimated the amount of flux 
emitted at IR wavelengths using the behaviour of SNe 2006aj, 
2009bb and 2010bh as a proxy. Note also that the bolometric 
LCs of SN 19941 and SN 2003dh are in filters BVR C I C and 
UBVR C I C respectively. Finally, usin g the same method that 
we applied to SN 2010bh in Section I6TT1 we have estimated 
the amount of flux emitted at UV wavelengths for SN 1998bw 
by using SN 2006aj as a template/proxy. 

As for SN 2006aj, we fit our analytical model to the bolo- 
metric light curves of the comparison SNe. We have used 
values of the photospheric velocities given in the literature for 
each event. The results of our fits are presented in Table |6] 
where the quoted errors are statistical only. 

When we compare the values of the nickel mass (MnO, 
ejected mass (M e j) and explosion energy (E^) determined 
from our model with those in the literature, it is seen that 
our model is able to recover the physical parameters of these 
events that are determined by more detailed spectral and pho- 
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FIG. 9.— The quasi-bolometric UBVR C UH light curves of XRF 100316D / 
SN 2010bh as well as three spectroscopically-confirmed GRB-SNe: 1998bw, 
2003dh & 20()6aj and two local type Ic SNe that were not accompanied by 
a GRB-trigger: Ic SN 19941 & broad-lined Ic SN 2009bb. All of the bolo- 
metric LCs have been corrected for foreground and host extinction. Note 
that the bolometric LC of SN 19941 and SN 2003dh are in filters BVR C I C and 
UBVR C I C respectively. 

tometric modeling. We note that on average, we tend to 
under-estimate the amount of ejected mass, which in turn 
under-estimates the explosion energy in these events, apart 
for SN 2006aj, where we over-estimate the explosion energy. 
There are several assumptions that go into our model, such 
as no asphericity and the central location of the nickel mass 
during the explosion (that does not mix with the expanding 
ejecta), as well as a constant value for the optical opacity. 
Even if the opacity is approximately constant during the pho- 
tospheric phase, it may be a slightly different value to that we 
are using here. 

Finally, the nickel mass determined from our fit of 
SN 2010bh (Section 5) reiterates the faint nature of 
SN 2010bh in relation to the other GRB-SNe. While the 
nickel mass produced by SN 2010bh is slightly greater than 
that of SN 19941, it is less than that seen in other GRB-SNe. 
Indeed, for all of the GRB-SNe that have been investigated 
through spectral and photometric modeling, SN 2010bh ap- 
pears to be the faintest (i.e. has created the least amount of 
nickel during the explosion). 

In conclusion, the energetics of GRB-SNe and Ic-BL SNe 
such as SN 2009bb and SN 2002ap are comparable, and 
are considerably higher than that of other Ibc SNe such as 
Ic SN 19941 and lb SN 1999ex. Additionally, the ejected 
masses of the GRB-SNe and Ic-BL SNe are commensurate 
and somwhat higher than some Ibc SNe (e.g. SN 19941 and 
SN 1999ex), though some Ibc SNe do ejected comparable 
masses (e.g. SN 2009jf and SN 1999dn). Finally, the nickel 
masses of all these events span a range of 0.07 — 0.52M E ,, 
with no clear distinction between GRB-SNe, Ic-BL SNe and 
Ibc SNe. 

6.3. Stretch and Luminosity Factors of the GRB-SNe 

In an attempt to put the luminosity and stretch factors of 
SN 2010bh (Table |5]l into context, we assembled from the lit- 
erature all of the photometric data available for most of the 
previously detected GRB-SNe. Note that we have not pur- 
sued every event but only included those where several mea- 
surements were made of the supernova "bump" that would 
allow us to constrain the properties of the associated SN. 
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TABLE 6 

Physical parameters of Ibc SNe and GRB-SNe 



SN 


type 


Vph (km s ) 


M Ni (Mr.) 


M eJ (M Q ) 


E k (10 52 erg) 


Wavelength Range (A) 


Ref. 


Note" 


19941 


Ic 


10,000 


0.07 ±0.01 


0.88 ±0.05 


0.09 ±0.01 


4 4()0 - 


- 8 000 


(1) 


This Work 


19941 


Ic 


10,000 


ss 0.07 


rsQ.9 


0.10±0.01 


4 400 - 


- 8 000 


(12) 


Literature 


GRB 980425 / 1998bw 


Ic-BL 


18,000 


0.43 ± 0.02 


7.04 ±0.57 


2.27±0.17 


3 000 — 


16 600 


(2) (3) (4) (5) 


This Work 


GRB 980425 / 1998bw 


Ic-BL 


18,000 


0.4-0.7 


8±2 


2-5 


3 000 — 


16 600 


(13), (14) (5) 


Literature 


GRB 980425 / 1998bw 


Ic-BL 


18,000 


0.52 ±0.02 


5.36±0.53 


1.73 ±0.17 


1 QOO 


16 600 


(2) (3) (4) (5) 


This Work 


1999dn 


lb 




RS 0.11 


4-6 


0.5 


^ ooo 


1 fx f\(Y\ 

io, ouu 


(18) 


Literature 


1999ex 


lb 




0.25 


0.9 


0.03 






(16) 


Literature 


2002ap 


Ic-BL 




ss 0.07 


2.5-5.0 


0.4- 1.0 


-J, uuu 


1 fi ^nn 
io, ouu 


(17) 


Literature 


GRB 030329 / 2003dh 


Ic-BL 


20,000 


0.39 ±0.03 


4.82 ±0.54 


1.92±0.21 


t, 000 

J , uuu - 


s nnn 

- o , uuu 


(6) 


This Work 


GRB 030329 / 2003dh 


Ic-BL 


20,000 


S3 0.4 


7±3 


3.5± 1.5 


^ 000 

-J , uuu 


s nnn 

o , uuu 


(6) 


Literature 


2003jd 


Ic-BL 




0.36 ±0.04 


3.0 ±0.5 


ks 0.7 


4,400- 


-8,000 


(19) 


Literature 


2004aw 


Ic 




0.2-0.3 


3-8 


«0.8 


3,000- 


16,600 


(19), (20) 


Literature 


XRF 060218 / 2006aj 


Ic-BL 


20,000 


0.14 ±0.01 


1.92±0.12 


0.76 ±0.05 


3,000- 


16,600 


(7), (8), (9), (10) 


This Work 


XRF 060218 /2006aj 


Ic-BL 


20,000 


a 0.2 


sa2 


«0.2 


3,000- 


16,600 


(15) 


Literature 


XRF 060218 / 2006aj 


Ic-BL 


20,000 


0.16 ±0.01 


1.63 ±0.05 


0.65 ±0.02 


1,900- 


16,600 


(7), (8), (9), (10) 


This Work 


2009bb 


Ic-BL 


15,000 


0.23 ±0.01 


3.86±0.13 


1.73 ±0.06 


3,000- 


16,600 


(11) 


This Work 


2009bb 


Ic-BL 


15,000 


0.22 ±0.06 


4.1 ±1.9 


1.8±0.7 


3,000- 


16,600 


(11) 


Literature 


2009jf 


lb 




0.17 ±0.03 


4-9 


0.3-0.8 


3,000- 


-8,000 


(21) 


Literature 


XRF 100316D/2010bh 


Ic-BL 


25,000 


0.10 ±0.01 


2.24 ±0.08 


1.39 ±0.06 


3,000- 


16,600 


This Work 




XRF 100316D/2010bh 


Ic-BL 


25,000 


0.12 ±0.01 


1.93 ±0.07 


1.20 ±0.05 


1,900- 


16,600 


This Work 





° Physical parame ters determined in: (1) this paper, or (2 ) taken from the l iteratu re. 
(1) Richmond et al. 1 1996), (2) Galama et al. 1 1998a), (3) McKenzie & Schaefei ( 1999), (4) Sollerman et al. (2002), (5) Patat et al. (2001), (6) Deng et al. 
(2005), (7) Campana et al. (2006), (8) Sollerman et al. (2006), (9) Cobb et al. (2006), (10) Ferrero et al. (2006), (11) Pignataetal. (2011), (12) Iwamoto et al. 
(1994), (13) Iwamoto et al. 1 1998), (14) Nakamura et al. 12001), (15) Mazzali et al. (2006), (16) Richardson et al. 12006) and references therein, (17) 
Mazzali et al. (2002), (18) Benetti et al. (2011), (19) Valenti et al. (2008), (20) Taubenberger et al. (2006), (21) Sahu et al. (2011) 



For every GRB-SN event we assume that light is coming 
from three sources: (1) the afterglow, (2) the supernova and 
(3) the host galaxy. First we remove the constant source of 
flux from the host galaxy, either by using magnitudes deter- 
mined by image subtraction by the authors (e.g. XRF 020903; 
Bersier et al. 2006), or by mathematically subtracting the flux 
from each measurement using observations made of the host 
galaxy long after the OT had faded below the brightness of 
the galaxy (e.g. GRB 090618; Cano et al. 201 1). 

Then, we use the early-time data for each event, and model 
the behaviour of the afterglow using single or broken power- 
laws: 

flux(t)=f X ((■=!— \ + ( — — ) ) (7) 

V V Tbreak J \ Tbreak J J 

where f is the flux zero-point, unique to each filter, t is the 
time since the burst and Tbreak is the time when the power- 
law changes from temporal index a,\ to a 2 - The afterglow 
parameters f , Tb rea k, OL\ and a 2 were allowed to vary in the 
fit for each event, and the results of our modeling are listed in 
Table □ 

In all cases the values we derived are similar to those seen 
in the literature. For example, for GRB 990712 we find that a 
single power-law fits both the V and /?-band data well with a 
decay constant of a = —0.95 ±0.02. In comparison, Bjorns- 
son et al. (2001) find that the V-band data can be fit with either 
a single or broken power-law, where the value in the former is 
ay — 0.82 ± 0.03. They find that the /?-band photometry can 
be fit with a single power-law with a slightly steeper decay 
index of ccr = —0.91 ±0.06. Both of these values are fully 
consistent with our own fit of the optical data. 

Additionally, for GRB 020405 Price et al. (2003) found 
from their /?-band photometry a decay constant of a — — 1 .41 



when they limited the range of the fit up to 10 days after the 
burst, while Masetti et al. (2003) find a = - 1 .54 ± 0.06 in all 
optical filters for the same time period. However, when Price 
et al. (2003) re-fit the data with a broken power-law, they find 
a better fit with parameters a x = -0.94 ± 0.25, a 2 = - 1 .93 ± 
0.25, and the time the LC breaks from <X\ to a 2 of Tbreak — 
1.67 ±0.52 days (x 2 /d.o.f= 35.9/28). Similarily, Masetti et 
al. (2003) find that at times > 20 days, a = -1.85 ±0.15. 
When we restrict our fit to 2 — 10 days, we find that a single 
power-law provides a suitable fit to the data with a = 1.72± 
0.08, j 2 /d.o.f = 6.16/4. This value for the temporal index is 
consistent with the value of a.2 found by Price et al. (2003) 
and the value of a found by Masetti et al. (2003) at late times. 

Once the afterglow parameters were determined, we then 
subtract Equation (|7]i from the OT (optical transient, e.g. af- 
terglow and supernova) flux light curve to obtain host- and 
afterglow-subtracted "supernova" light curves (e.g. Cano et 
al. 2011). The resultant l ight curves are displayed in Fig- 
ured)] Note that in Figure[K)]for XRF 060218 / SN 2006aj, 
we have included an additional power-law component of a = 
—0.2, which was fixed in the fit, and is similar to that deter- 
mined by Ferrero et al. (2006) of a = -0.20 ±0.10. 

For each of the host- and afterglow-subtracted "supernova" 
light curves, we determined the stretch and luminosity factor 
for each GRB-SN in relation to SN 1998bw in various fil- 
ters, repeating the same process we followed for SN 2010bh 
(Section 3.5). The results of our fit are listed in Table [8] with 
the /?-band stretch and luminosity factors displayed in Figure 
[TT1 and the host- and afterglow-subtracted "supernova" light 
curves displayed in Figure [TOl 

A previous study by Zeh et al. (2004), which was ex- 
tended by Ferrero et al. (2006) (F06 hereon), also determined 
the stretch and luminosity factors of GRB-SN in relation to 
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FIG. 10. — Mosaic of GRB-SNe. For each event the host has been subtracted, either by image subtraction by the original authors or by mathematically 
subtracting the flux due to the host. The OT (afterglow and supernova) light curves were then modeled, with the parameters of the afterglow determined and 
then subtracted out to create host- and afterglow-subtracted "supernova" light curves. Plotted for each event (dashed line) is how SN 1998bw would appear if it 
occurred at the redshift of the event. The solid-line is the light curve of SN 1998bw corrected by the luminosity and stretch factors listed in Table [8] For XRF 
060218 / SN 2006aj, we have included an additional power-law component to the SN with index a = —0.2, similar to that found by Ferrero et al. (2006). 



XRF 100316D / SN 2010bh and the Nature of Gamma Ray Burst Supernovae 



15 



2 - 



1 - 
0.8 - 
0.6 " 



0.1 - 




0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



s 

FIG. 11. — Luminosity factor (k) versus stretch factor (s) in the R band 
for all of the GRB-SNe in Tables \5\ and [8] We find a correlation/Pearson 
coefficient for the entire sample of c = —0.05, with a 95% confidence interval 
of (—0.46,0.38) for a significance of p > 0.05. As the critical value for the 
Pearson correlation coefficient for N = 22, d.o.f. = 20 and p > 0.05 is 0.492, 
we can conclude there is not a statistically significant correlation. 

SN 1998bwE3by modeling the SN "bumps" seen in the GRB 
light curves. For many events where we have determined the 
stretch and luminosity factors from the afterglow- and host- 
subtracted "supernova" light curves, and F06 from the super- 
nova "bumps", we find similar results: e.g. GRB 050525A 
/ SN 2005nc: this paper: k = 0.69 ±0.03, s = 0.83 ±0.03; 
FOfQ k = 0.66^, s = 0.77 ±0.04; and XRF 060218 / 
SN 2006aj (V band): k = 0.72 ±0.02, s = 0.65 ±0.02; F06: 
k = 0.724 ± 0.007, s = 0.682 ± 0.005, with similar results for 
the other filters. 

However, for some events we find different properties of 
the GRB-SN: for GRB 01 1 121 / SN 2001ke, we find similar 
stretch factor to F06 (this paper: s = 0.84 ±0.01; F06: s = 
0.80 ± 0.02), however we find SN 2001ke to be much brighter 
once correcting for host-extinction: (this paper: k = 2.11 ± 
0.03; F06: k = 0.88^;^). Additionally, for GRB 020405 
we find a similar luminosity factor as F06: (this paper: k = 
0.82 ± 0. 14; F06: k = 0.90+g ;J^), however we find that the SN 
evolves more quickly: (this paper: s — 0.62 ±0.03; F06: s — 
0.97 ± 0.07). However, GRB 020405 aside, for every other 
GRB-SN event we find the same value for the stretch factor 
as F06 within our respective errorbars. 

We note that when using this method to determine the lu- 
minosity factor relative to SN 1998bw, the SN brightness de- 
pends on how accurately we have modeled the afterglow. For 
instance, if we over-estimate the relative contribution of the 
afterglow, then we will correspondingly under-estimate the 
SN brightness. 

We find that the stretch factor of each GRB-SN event is 
approximately the same in each filter (e.g. GRB 9907 12, GRB 
031203, XRF 060218, GRB 060729 & GRB 090618), though 
this trend applies less in the bluer filters. For example, for 
XRF 060218 / SN 2006aj, we find stretch factors in the R 
and / filters of s = 0.68 ± 0.02, while in V and B we find 
slightly smaller values of s = 0.65 ± 0.02 and s = 0.60 ± 0.02 
respectively. 

We have investigated whether there might be a correlation 
between the k and s factors of the GRB-SNe, perhaps some- 

23 only in the R band 

24 in the R band and corrected for host-extinction 



thing akin to that established for type la SNe (i.e. the possibil- 
ity of using GRB-SNe as "standard candles"). We have used 
all of the fi-band data, of which there are 18 GRB-SNe. For 
GRB/XRFs 02 1 2 1 1 , 040924, 0803 1 9B and 09 1 1 27, for which 
no R band data exists, we have assumed that the stretch and 
luminosity factors in the V band (GRB 02 1 2 1 1 ) and the / band 
(GRBs 040924, 0803 19B and 091 127) are approximately the 
same as in the R band. As already noted for events such as 
XRF 060218 and GRB 090618, the stretch and luminosity 
factors in filters VR C I C are approximately the same and differ 
by only a few percent. The inclusion of these events brings 
the total sample size up to N = 22, and for the entire sample 
we find a correlation/Pearson coefficient of c = —0.05 and a 
95% confidence interval of (—0.46,0.38) for a significance of 
p > 0.05. As the critical value for the Pearson correlation co- 
efficient for N = 22, d.o.f. = 20 and p > 0.05 is 0.492, we 
can conclude that a statistically significant correlation does 
not exist. 

Stanek et al (2005) noted that a relation may exist between 
light curve shape and luminosity, which the results of F06 
could neither confirm nor refute (total sample of N = 13). 
F06 noted that there might be a trend of rising k with rising s, 
however the addition of 9 more GRB-SNe to the total sample 
shows that it is highly unlikely that such a correlation exists. 
When this result is taken in tandem with the results of Drout 
et al. (2011), who found no statistically-significant correla- 
tion between peak luminosity and light-curve decay shapes 
for a sample of Ibc SNe, indicates that Ibc SNe (including 
GRB-SNe and Ic-BL SNe) cannot be used as standardizable 
candles. 

7. SUMMARY 

We have presented optical and near infrared photometry of 
XRF 100316D / SN 2010bh obtained on the Faulkes Tele- 
scope South, Gemini South and HST, with the data span- 
ning from t — 1 = 0.5 —47.3 days. It was shown that the 
optical light curves of SN 2010bh evolve more quickly than 
the archetype GRB-SN 1998bw, and at a similar rate to 
SN 2006aj, which was associated with XRF 060218, and non- 
GRB associated type Ic SN 19941. In terms of peak lumi- 
nosity, SN 2010bh is the faintest SN yet associated (either 
spectroscopically or photometrically) with a long-duration 
GRB or XRF, and has a peak, V-band, absolute magnitude 
of M v = -18.62 ±0.08. 

SN 2010bh appears to be redder than GRB-SNe 1998bw 
and 2006aj, where the colour curves of SN 2010bh are seen 
to be redder early on, though at late times the V —R and R — i 
colour curves matched that of SN 19941. The red nature of 
SN 2010bh is also demonstrated in Figure [8] where it is seen 
that more of the bolometric flux is emitted at infrared wave- 
lengths (and less at optical wavelengths) than in the broad- 
lined Ic SN 2009bb. 

We also gave evidence of the detection of light coming from 
the shock break-out at t — t G = 0.598 days. The brightness 
of the B-band light curve at this epoch, as well as the shape 
of the optical SED, which with a spectral power-law index 
of j3 = +0.94 ± 0.05 (which is harder than is expected for 
synchrotron radiation), implies that the source of light at this 
epoch is not synchrotron in origin and is likely coming from 
the shock-heated, expanding stellar envelope. 

We then applied a simple physical model to the bolometric 
light curve of SN 2010bh. When we include all photometry 
in the optical and infrared regime (in the 3, 000 A — 16, 600 A 
wavelength range) we find physical parameters of Mm; = 
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TABLE 7 

Afterglow parameters of the GRB-SNe 



GRB/XRF/SN 


Filter 


Oft 


«2 


Tbreak (days) 


*7d.o.f 


011121/2001ke 


R c 


-1.70±0.01 






3.56/4 


020405 


Rc 


-1.72±0.16 






6.16/4 


020903 


Rc 


-1.87 ±0.05 






1.86/2 


041006 


Rc 


-0.47 ±0.35 


-1.47 ±0.39 


0.30 ±0.78 


3.24/3 


050525 A/2005nc 


R c 


-0.89±0.22 


-2.13±0.30 


0.43 ±0.35 


4.68/3 


060729 


u * 


-0.01 ±0.03 


-1.65 ±0.05 


0.75 ±0.08 


37.99/29 


090618 


Rc, i f 


-0.65 ±0.07 


-1.57 ±0.07 


0.50 ±0.11 


151.42/113 



* Cano et al. (2011) determined the afterglow parameters using the {/-band data due to the 
lack of early-time data in R and /. 

' Cano et al. (201 1) fit all of the optical data simultaneously. 



TABLE 8 

Luminosity (it) and Stretch (.«) Factors of the GRB-SNe with respect to 

SN 1998BW 



GRB/XRF/SN 


Filter 


s 


k 


Ref (s & k) 


Ref (Photometry) 


970228 


R 


1.45 ±0.95 


0.40 ±0.29* 


(1) 




980425 


BVRi 


1.00 


1.00 




(4) 


990712 


V 


0.83 ±0.44 


0.65 ±0.46 


This paper 


(5), (6), (7), (8) 


990712 


R 


0.90 ±0.07 


1.45 ±0.20 


This paper 


(5), (6), (7), (8) 


991208 


R 


1.12±0.20 


2.62±1.10 


(1) 




000911 


R 


1.40 ±0.32 


0.85 ±0.44 


(1) 




011121/2001ke 


R 


0.84 ±0.01 


2.11 ±0.03 


This paper 


(9), (10) 


020405 


R 


0.62 ±0.03 


0.82±0.14 


This paper 


(11), (12) 


020903 


R 


0.85 ±0.27 


0.71 ±0.40 


This paper 


(13) 


02121 l/20021t 


V 


0.98 ±0.26 


0.40 ±0.19* 


(1) 




030329/2003dh 


R 


0.81 ±0.02 


0.90 ±0.04 


This paper 


(14), (15) 


031203/20031w 


R 


1.08 ±0.02 


1.51 ±0.03 


This paper 


(16), (17), (18) 


031203/20031w 


I 


1.07 ±0.04 


1.29 ±0.04 


This paper 


(16), (17), (18) 


040924 


1,1 


1.00 (fixed) 


«0.25 


(2) 




041006 


R 


1.47 ±0.04 


1.16 ±0.06 


This paper 


(19) 


050525A/2005nc 


R 


0.83 ±0.03 


0.69 ±0.03 


This paper 


(20) 


050824 


R 


0.52±0.14 


1.20 ±0.69 


(23), (24) 




06021 8/2006af* 


B 


0.60 ±0.01 


0.67 ±0.02 


This paper 


(1) 


06021 8/2006aj** 


V 


0.65 ±0.01 


0.72±0.01 


This paper 


(1) 


060218/2006af* 


R 


0.68 ±0.01 


0.72 ±0.01 


This paper 


(1) 


06021 8/2006af* 


1 


0.68 ±0.01 


0.76 ±0.01 


This paper 


(1) 


060729 


R 


0.89 ±0.04 


0.84 ±0.05 


This paper 


(21) 


060729 


I 


0.94 ±0.02 


1.03 ±0.02 


This paper 


(21) 


0803 19B 


I 


0.89 ±0.10 


2.30 ±0.90 


(3), (25) 




090618 


R 


1.09 ±0.04 


0.70 ±0.05 


This paper 


(21) 


090618 


I 


1.06 ±0.03 


0.94 ±0.03 


This paper 


(21) 


091127 


I 


0.88 ±0.01 


0.89 ±0.01 


This paper 


(22) 


101225A 


R,i,Z 


1.25 ±0.15 


0.08 ±0.03 


(26) 





* Not corrected for host-extinction. 

** Fit excluding additional decaying power-law component (see text). 

(DIFerrero et al.ld20 06). (2) So derberg et"aTlll2006l). (3)ITanvir et alj j20T0h. (4)IGalama et al I 
<1998al). (5) ISahu et alj 12000). (6) Christensen et al. (2004), (7) Biornsson et al.1 1200111 . 
(8) Hiorthetal. (2000), (9) Garnavich et al. (2003), (10) Kiipcii Yoldas et al. (2007), (11) 
IBersier et al.l J2003I). (12) IPrice et alj (120031), (13) IBersier et alH 2006). (14) Mathes mTetaT] 
120031). (15)beng et alj (2005), (16) Malesani et al. (2004), (17) Mazzali et al. l200 3), (18) 
Margutti et al. (2007), (19) Stanek et al. (2005), (20) Delia Valle et al. (2006), (21) Cano et al. 
<2TjTl|). (22) ICobb et al l <20ToT), (23) Sollerman et al. (2007), (24) Kann et al. 1 20lQ), (25) 
Bloom et al. (2009), (26) Thone et al. l20n|) 
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O.1O±O.O1M , M ej = 2.24±O.O8M , E k = 1.39 ± 0.06 x 
10 52 erg. The faint nature of SN 2010bh becomes again ap- 
parent when the nickel mass synthesized during the explo- 
sion is compared with other GRB-SNe such as SN 1998bw, 
where it is believed that Mnj « O.4M was created in the ex- 
plosion. Indeed SN 2010bh synthesized only a marginally 
larger amount of nickel than local type Ic SN 19941 (Mm; k 
O.O7M ). The investigation of SN 2010bh in relation to the 
general population of type Ic SNe has once again shown that 
type Ic SNe are a very heterogeneous class of supernovae, 
spanning a wide range of luminosities and ejected masses. 

Finally, we assembled from the literature all of the avail- 
able photometry of previously detected GRB-SNe. For each 
of these events we assumed that light is coming from three 
sources: from the host galaxy, the afterglow and the super- 
nova. First we removed the constant flux due to the host 
galaxy, then we modeled the optical afterglows, and then 
subtracted the light due to the afterglow, which resulted in 
host- and afterglow-subtracted "supernova" light curves (Fig- 
ure [Toll. We then created synthetic light curves of SN 1998bw 
as it would appear if it occurred at the redshift of each given 
GRB-SNe. We then compared the brightness and shape of the 
GRB-SNe with that of the shifted LCs of SN 1998bw, and 
derived stretch (s) and luminosity (k) factors for each GRB- 



SNe. The results of our method, when compared to previous 
studies performed by Ferrero et al. (2006) showed similar val- 
ues for the stretch and luminosity factors, as well as including 
more events than the previous studies. When we checked for 
the possibility of a correlation between k and s in the /?-band 
data, we find an insignificant correlation/Pearson coefficient 
of c = —0.05, suggesting that it is highly unlikely that such a 
relation exists. 
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